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Abstract 



^ ' A Roy equation analysis of the available tttt phase shift data is 

5-H ' 



performed with the I = S- wave scattering length a{] in the range 
predicted by the one-loop standard chiral perturbation theory. A suit- 
able dispersive framework is developed to extract the chiral coupling 
constants l\, I2 and yields l\= — 1.70±0.15 and 5.0. We remark on 
the implications of this determination to (combinations of) threshold 
parameter predictions of the three lowest partial waves. 
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1 Introduction 



Chiral perturbation theory ]l|, |2|, |3| provides the low energy effective theory of 
the standard model that describes interactions involving hadronic degrees of 
freedom and exploits the near masslessness of the u, d (and s) quarks and the 
observation that the pions, kaons and the rj could be viewed as the Goldstone 
bosons of the spontaneously broken axial symmetry of massless QCD. It is 
a non-renormalizable theory and involves additional coupling constants that 
have to be introduced at each order in the derivative or momentum expansion. 
[From here on we confine our attention to the SU{2) flavor subgroup.] At 
leading order, 0(p 2 ), we have the pion decay constant, F n ~ 93 MeV in 
addition to the mass of the pion, m n = 139.6 MeV, henceforth set equal to 1. 
As a result, one has for what is arguably the simplest purely hadronic process 
of 7T7T scattering a prediction for a key threshold parameter, the / = 0, S- wave 
scattering length a[| = 7 / (32irF%) ~ 0.16 [Q. There are ten more coupling 
constants at the next to leading order, 0(p 4 ); four of them, li, l 2 , I3 and U 
enter the nir scattering amplitude ffl. As a result, at this order a° (and Oq) 
are predicted in terms of these as well. One of the cornerstones of standard 
chiral perturbation theory at 0(p A ) is a relatively definite prediction of a[] in 
the range 0.19-0.21. 

The coupling constants l\ and I2 have been fixed in the past from ex- 
perimental values available in the literature || for the D- wave scattering 
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lengths: 



^ = 17 ± 3 • 1(T 4 , a\ = 1.3 ±3- 1(T 4 . 



The one-loop expressions for these are 

4 = ^-r(Zi + 4r 2 -— ), a\ = 1 —-(l 1 + l 2 -—) (l.i) 

2 1440vr 3 Fi V1 2 8 ; ' 2 14407r 3 Fi V1 2 40 ; V ; 



and yield 

Zi = —2.3 ±3.7, f 2 = 6.0 ±1.3. 
These have also been determined from analysis of decays M which yield 

h = -0.7 ± 0.9, l 2 = 6.3 ± 0.5 
and more recently by estimating higher order corrections to these decays [0 

f x = -1.7 ±1.0, r 2 = 6.1±0.5. 

l\ and I2 are coupling constants consistent with the presence of resonances. In 
particular the p resonance may make a significant contribution [see Appendix 
C in Ref. ||] and has also been discussed more recently ||. Furthermore, 
generalized vector meson dominance leads to numerical values for these 
consistent with the numbers above. Tensor resonances also have been found 



to make non-trivial contributions 1 10 



The constants I3 has been estimated from the analysis of SU(3) mass 
relations which yields 0: 

l 3 = 2.9 ± 2.4. 
3 



The variation of a[J is essentially equivalent to the variation of I3. While 
here I3 would have to be —70 in order to accommodate a[j = 0.26, there 
is an extended framework which re-orders the chiral power counting in or- 
der to accommodate such large values of aj] modifying even the tree-level 
prediction Jll]]. Here we confine ourselves to the more predictive standard 
chiral perturbation theory whose stringent predictions will come under ex- 
perimental scrutiny Thus we note that in our final analysis we cannot 
claim an independent determination of I3 via the Roy equation analysis of 
this work since a[j is varied anyway in the range predicted by standard chiral 
perturbation theory. 

The constant U enters the one-loop expression of the relatively accurately 
determined value of the "universal curve" combination 2a[] — 5<2q Jl3j and is 



also related to the independent estimates of the scalar charge radius of the 
pion. An SU (3) analysis that has been performed for the ratio of the kaon 
and pion decay constants Fk/F w also provides a measure of this constant 
Z4 ~ 4.6 [0. In the present analysis we treat a[] as the only free parameter 
to the fit to the data and is produced as an output corresponding to the 
optimal solution of our data fitting. In particular, the values we find remain 
consistent with the universal- curve band. Thus we have a determination of 
the constant 1^. However, we also perform constrained fits with computed 
from certain universal curve relations that fix Z 4 a priori. The influence on 
the actual numerical fits is found to be minimal reflecting the weakness of 



the 1 = 2 channel and influences the determination of 1% and I2 minimally 
due to reasons we discuss in subsequent sections. 

On the other hand irn scattering has been studied in great detail in ax- 
iomatic field theory fll5"l . (Fixed-t) dispersion relations with two subtractions, 
a number dictated by the Froissart bound, have been rigorously established 
in the axiomatic framework. The properties of crossing and analyticity have 
been exploited in order to establish the Roy equations |16|, [L7]], a system of 
integral equations for the partial waves. The Roy equations have been the 



basis of analysis of phase shift data |L8 and a knowledge of the threshold 



parameters involved in mr scattering has been obtained. Best fits to Roy 
equation analysis of data are obtained with a[] = 0.26 ± 0.05 [ l9[ . Note that 



the D- wave scattering lengths cited earlier have also been extracted from 
Roy equation analysis. The properties of analyticity, unitarity and crossing 
and positivity of absorptive parts have also been shown to produce non-trivial 
constraints on the magnitudes of a certain combination of l\ and I2 f20] . 

Here we report on a direct determination of the chiral coupling constants 
from the existing phase shift data itself by performing a Roy equation fit to 
it when a[j is confined to the range predicted by chiral perturbation theory. 
The chiral amplitude is rewritten in terms of a dispersive representation with 
a certain number of effective subtractions consistent with 0(p 4 ) accuracy, 
where the subtraction constants are expressed in terms of the chiral coupling 
constants. The fixed-t dispersion relations of axiomatic field theory are also 



rewritten in a form whereby a direct comparison can be made with the chiral 
dispersive representation, while the effective subtraction constants are now 
computed in terms of physical partial waves produced by the Roy equation 
fit, the input value of a[j and the resulting value of generated by the fit. In 
most of our treatment we limit ourselves to a certain approximation where 
we account for the absorptive parts of I > 2 states only through the "driving 
terms" of the Roy equations for the S- and P- waves. 

Furthermore, we also perform an analysis of certain threshold parame- 
ters computed from the Roy equation fits which reveals the magnitude of 
0(p 6 ) corrections their one- loop predictions are expected to suffer from. The 
dispersive framework can be extended to meet the needs of two-loop chi- 
ral perturbation theory and used to pin down the associated coupling con- 
stants pi], |22|. The work reported here summarizes the first stage of our 
computations and is presently being extended to meet the needs of the two- 
loop computation of P2f . 
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2 7T7T Scattering to 0(p A ) in chiral perturba- 
tion theory and the Roy equation solutions 

The notation and formalism that we adopt in this discussion follows that of 



Ref. fla] . Consider titt scattering: 

^ a (Pa)+7T b ( Pb )^TT C (p c )+7r d ( Pd ), 

where all the pions have the same mass. The Mandelstam variables s, t and 
u are defined as 

S=(Pa+Pb) 2 , t=(Pa-Pc) 2 , t=(Pa-Pd) 2 , S + t + U = 4. (2.1) 

The scattering amplitude F(a, b — *> c,d) (our normalization of the amplitude 
differs from that of Ref. [18| by 32ir and is consistent with that of Ref. P, |23fl) 
can then be written as 

F(a, b — > c,d) = 5 ab 5 cd A(s, t, u) + 5 ac 5bdA(t, s, u) + 5 ad 5 bc A(u, t, s). 

From A(s,t,u) we construct the three s-channel isospin amplitudes: 

T°(s,t,u) = 3A(s,t,u) + A(t,s,u)+A(u,t,s), 

T^s^u) = A(t,s,u)-A(u,t,s), (2.2) 

T 2 (s,t,u) = A(t,s,u) + A(u,t,s). 

We introduce the partial wave expansion: 

T\s,t,u) = 327^(2/ + mC—^fRs), (2.3) 
i=o s 4 



fi(s) = ft{s) = 0, 1 odd, fl(s) = 0, 1 even. 

The unitarity condition for the partial wave amplitudes f{(s) is: 

toff {*) = p(s)|//( S )| 2 + 1 " ( f ( , a))a , (2.4) 

4p(s) 



where p(s) = J(s — A)/s and 77/ (s) is the elasticity at a given energy a/s. 
We also introduce the threshold expansion: 

Re/fr) = (^) ' (af + 6f (^) +...), (2.5) 

where the af are the scattering lengths and the 6f are the effective ranges, 
namely the leading threshold parameters. 

Chiral perturbation theory at next to leading order gives an explicit rep- 
resentation for the function A(s,t,u) at 0(p 4 ) [0]: 

A(s,t,u) = A {2) {s,t,u) + A {4) {s,t,u) +0{p 6 ), (2.6) 



with 



s - 1 



A^(s,t,u) = ^(3(s 2 -l)J(s) 

+ [t(t -u)-2t + Au - 2)J(t) + (t++u] 
+ ^r^Wi - 4/3) (a - 2) 2 + (l 2 - 5/6) [s 2 + (t - it) 2 ] 
+12s(l 4 - 1) - 3(l 3 + AU - 5)} 
andJ(z) = l — f 1 dx lnfl - x(l - x)z], toJ(z) = ^-Q(z - A). 

I67H JO ' 167T 



Note also that at 0(p 4 ) the imaginary parts of the partial waves above thresh- 
old computed (s > 4) from the amplitude above is: 

l m/ 0( s ) = ( 2s -if 

JoK ' 1024vr 2 F/ V ; 

l m/ i( s ) = p( ' } A s -A) 2 

l m /2( s ) = A s-2) 2 (2.7) 

Im//(s) = 0, l>2. 



[Note that the chiral power counting enforces the property that the absorptive 
parts of the D- and higher waves arise only at 0(p 8 ).] Furthermore these 
verify the property of perturbative unitarity, viz., when the 0(p 2 ) predictions 
for the threshold parameters Oq = 7 /(32nF n 2 ), a 2 , = — 1 / (I6ir F n 2 ) , b^ = 
l/(47rF 7r 2 ), bl = -l/(87rF 7r 2 ) and a\ = l/(24vrF 7r 2 ) are inserted into the 
pertinent form of the perturbative unitarity relations: 

Im/ 7 ( S ) = p(s)(a I + ^( s -4)/4) 2 , 1 = 0,2 
lmfl(s) = p(s)(a\(s-A)/A) 2 . 

In order to carry out the comparison between the chiral expansion and the 
physical scattering data, we first recall that up to 0(p 6 ), it is possible to 
decompose A(s, t, u) into a sum of three functions of single variables as fol- 
lows 121: 
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A(s,t,u) = 32tt 



\w (s) + |( a - + ^(a - t)W l {u) 



+\ (w 2 (t) + W 2 (u) - ^W 2 (s) 



One convenient decomposition of the chiral one-loop amplitude is: 



(2.8) 



W (s) 



32n 
4 



s- 1 



+ — -( S -1/2) 2 J( S ) 



1 



967r 2 Fi 



^2(Zi - 4/3) (s - 2) 2 + 4/3(Z 2 - 5/6) (s - If 



+12s(k- 1) - 3(Z 3 + 4Z" 4 - 5)) 



(2.9) 



W 2 (s) = 



5767rFi 



( S -4)J( S ), 



16tt 



1 ^-2) 2 J( S ) + i ^(Z 2 -5/6)( S -2) 2 



4Fi 



(2.10) 
,(2.11) 



where we note that this decomposition is not unique, with ambiguities in the 
real part only. We observe that the imaginary parts of these functions verify 
the relation: 

ImW^x) = Iraft(x), I = 0, 2 
ImW^x) = lmfl(x)/(x - 4), 



which may be used to demonstrate the following dispersion relations: 
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W (s) 



-1 + 72h + 48/2 - 27k - 108f" 4 - 864tt 2 F^ 2 
9216tt 3 F/ 
59 - 144/'! - 96/2 + 2I6/4 + l728ir 2 Fj 
+ 18432^7 S (2 - 12) 

-797 + 360/i + 240/ 2 2 s 3 Z" 00 rfx 



s +-/ -v vW >), 



184320?r 3 F 7r 4 tt 7 4 x 3 (x - s) 

WAs) = + — / —. 4r, ^fl(x), (2.13) 

KJ 13824tt 3 F/ ti h x 2 {x-A){x-s) JlKh v ; 

6/2 - 5 23 - 24/2 60Z 2 - 77 2 

^ 2(S) " 11527T3F/ + 4608vr3F/ S + 46080^/' {2M > 

+ - f-^— Im/ 2 (x). 

7T J 4 X 3 (x — S) 

We now reconstruct A(s, t, u) from this dispersive representation for the Ws 
to obtain: 

s - 1 -540 + 480/"i + 960/" 2 - I8O/3 - 720/" 4 

A(s,t,u) = tt H o — 3 

V ' ' ; F^ 2 57607r 2 F^ 4 

110 + 480/"i - 720/4 _ 163 - 120/i 2 

5760vr 2 F u 4 S ~ 5760vr 2 F 7r 4 S 

460 - 480/~ 2 20u(s - t) + 20i(s - u) 

+ r r-, nr . 9p 4 r "I" 



5760vr 2 F 7r 4 v ; 5760vr 2 F^ 4 
154-120| 

5760vr 2 F 7r 4 V ; V ; 

1 s 3 f 00 dx 



+32vr 11- r (lm/ °(x) - Im/ 2 (: 

7T J4 X^(x — SJ V 

+ 2^ " U)t nT x 2 (x-1)(x-4) Im/ll(x) 

+§(*-*)- r 2 , d w zv W! 1 ^) 

2 7T J4 X Z [X — U){X — A) 

l/t 3 f da; j ~2 / \ ^ 3 Z^ 00 dx 2 , N 

+ o _ / "a? IT Im /o x + — / -§7 W 

2 1iri4 r(a;-f) 7r J4 x d (x — u) 
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This is seen to be the sum of a polynomial of second degree in s, t and 
u and a dispersive piece. The problem associated with the non-uniqueness 
of the real part of the decomposition into the W's is eliminated by setting 
u = 4 — s — t upon which we obtain a second degree polynomial in s and t: 



The Roy equation fit allows us to obtain a representation only for the 
S- and P- wave absorptive parts, [with some effects of higher angular mo- 
mentum states absorbed into the driving terms (see Appendix A)]. Thus, 
a determination of the physical S- and P- wave absorptive parts, allows us 
to construct a set of crossing symmetric amplitudes |17|, |25[ from which we 
extract a representation for A(s, t,u) (see Appendix A for details): 




(2.16) 




(2.17) 
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+ 3 2 {S ~ t] T r x\x -%{x - A) lmfl {X) 

2 \n Ja x 6 {x — t) it Ji x 6 {x — u) 



where 








= «0 


4 r°° dx j 
7i h x{x - 4) m/o 




1 

7T . 


r^lmfUx) 1 = 
U x 6 




3 

7T . 


L x^x-4) lmf ^ X) 


«0 


= /3? 


= Pi = o. 



1mf£(x) + - ^lmfi(x) 7 = 0,2 

7T J4 X 1 

(2.18) 



We are now in a position to compare the two representations for A(s,t,u) 
namely the chiral representation eq.( |2.15| ) and the axiomatic representation 
eq. (|2.17p . These are formally equivalent, with the dispersive integrals in the 
former described by chiral absorptive parts whereas in the latter by the phys- 
ical S- and P-wave absorptive parts. For the chiral expansion to reproduce 
low energy physics accurately we now require the effective subtraction con- 
stants to match. Once more setting u = 4 — s — t yields the polynomial piece 



of the representation eq.( 2.17 ): 



P = -1287r(/?l+ 7 §)(*-£-f) + + ll-W\)s 2 

+My - -og + 8[3{ - 16^)3 + 16n(al + 16 T g). (2.19) 

A straightforward comparison of eq. ( |2.16| ) and eq. fl2.19| ) yields explicit 
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expressions for li, I2, I3 and Z4. In particular we have for l\ and I2' 

2Q 

l 2 = 24 7 r 2 F/( 4807r2 ^ 4 + 32 7 r(/3 1 1 + 7o 2 )). (2.21) 

The actual numerical values we find for li and Z2 are reported in a subsequent 
section. These have an interesting dependence on the actual physical phase 
shifts: one observes that in eq.( [2.20 ), the presence of 7°. As a result we 



can anticipate Zi to be influenced by the input for a{j. In contrast, I2 has no 
dependence on 7$ and depends almost totally on the P- wave contribution 
via fll, as a result of the weakness of the 1 = 2 channel which renders 7q 
negligible in comparison with (3\ (and 7q). Since the P- wave happens to 
be the best determined experimental quantity, even the Roy equation fits to 
it are not strongly influenced by the input value of a{j. Thus we expect a 
determination of I2 in this manner to be very stable. 

The values we find for Z3 and I4 from the procedure above are in rough 
agreement with the estimates provided in the introduction, but suffer from 
the fact that here the comparison also involves the 0(p 2 ); they are determined 
only after large cancelations occur and thus a determination of these are not 
expected to be very reliable. As noted earlier the determination in this 
manner of T 3 cannot be considered independent of the input a[j and so we do 
not report it here. U may be determined from the dispersive formulas here. 
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3 Implications to Some Threshold Parame- 



Specific combinations of threshold parameters appear on the left hand side 
when Roy equations for some of the lowest partial waves (and for instance 
their energy derivatives) are evaluated at threshold. On the right hand sides 
one has energy integrals over partial wave absorptive parts. These serve as 
sum rules for such combinations. Indeed expressions for these sum rules have 
been derived even before the Roy equation program (see e.g., ||26|| ). While 
there are several inequivalent methods of obtaining such sum rules, it has 
been shown that as long as one is confined to the absorptive parts of the S- 
and P- waves alone, each of these methods would yield identical results for 
the right hand sides (for a recent discussion, see f27|). 

We consider the Roy equations eq. ( |A.3| ) in the following limits: 



These may then be rearranged to yield the specific combinations of threshold 
parameters satisfying: 



ters 




(3.1) 





(3.2) 
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-9 (a; - 4)Tmft(x) + 5(x - 4)Im/ 2 (:r)} 
+ 12 lim -^-Re (£(s, A), 
Qbl + (2a° - 5a 2 ) = - 

IT J 4 (X(X - 4)) 2 

[2(x - 4)Im/ °(x) + 9(x - A)lmfl(x) (3.3) 
+ (7x - 4)lmf*(x) -3(x - 2) y/x(x - 4)(a 2 ) 2 } 
+24 lim -^-Re d2(s, A), 
18a} - (2a°-5a 2 ) U> ' 



186 



7T J4 (x(x — 4)) 2 

{-2(x - 4)Im/°(s) + 9(3x - A)lmfl(x) 

+ 5(x - 4)Im/ 2 (*)} + 18 s lim ^jf|ffi (3.4) 
i _ 16 r A dx 

11 ~ ~ J 4 {x{x-A)y 

{-2(x - 4) 3 Im/ (x) + 9(3x 3 - 12x 2 + A8x - 64)Im/ 1 1 (a;) 

+ 5(,-« W } + 72taA(^i|^)). (3 . 5) 

From the following limits for the Roy equations: 

, lm _MW , S „J« p.6) 

*-4+ (( s _4)/4) 2 ' s-4+ (( s - 4)/4) 2 V ; 

we find expressions of sum rules for the D-wave scattering lengths. 
16 r A dx 



a 2 = 



4577 J 4 x 3 (x — 4) 

{(x - 4)Im/ °(x) + 9(x + 4)Iin/ 1 1 (:r) + 5(x - 4)Iin/ 2 (+)} (3.7) 
Re d° 2 (s,A) 



+ lim 



s-4+ (( s _4)/4) 2 ' 

16 



907T J 4 x 3 (x - 4) 

{2(x - 4)Im/ °(x) - 9{x + 4)Im/ 1 1 (x) + (x - 4)Im/ 2 (s)} (3.8) 

+ hm - 2 \ ' / 

s^4+ (( s - 4)/4) 2 

The one-loop expressions for the combinations of interest are 

36 o _ (2fl o _ 5fl 2 } = __l__(2r i + 3r 2 - |), (3.9) 

66 2 + (2a° - 5a 2 ) = ^^(h + 3T a + jj), (3.10) 

1 55 
18a} - (2a° - 5a 2 ) = (3.11) 

1 - - 97 

186} = a I + ) (3.12) 

1 16vr 3 F4 V2 1 120 ; V ; 

respectively. Those for the D-wave scattering lengths are given in eq. ( |1 . 1| ) . 
Since the Roy equations are a consequence of dispersion relations with two 
subtractions, we note that in all the above the leading 0(p 2 ) contribution 
cancels exactly. Furthermore, the constants I3 and U are also absent since 
they accompany only constant and linear pieces in s, t and u in the 0(p A ) 
scattering amplitude. It would therefore seem that any two of the six combi- 
nations above is suitable for determination of l\ and I2 since they enjoy the 
same status as the D-wave scattering lengths (see eq. (|1.1|)). 

A careful consideration of these reveals some interesting characteristics. 
Note for instance that the one-loop expressions for 18a} — (2a[j — 5a 2 ,) and 
for 186} depends only on the combination h — h- If we now evaluate their 
numerical values for each of our Roy equation phase shift representations by 
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inserting them into the sum rules in order to compute the coupling constants 
of interest we would not have a meaningful result since each of the combina- 
tions above represents parallel straight lines in the Zi, Z 2 plane! If we consider 
the further combination 18a\ — (2a{] — 5<2q) — 18&J, its one-loop expression 
is Sl/l^On^F^. If we saturate the sum rule for this combination with 
A = 00 and the chiral absorptive parts eq.( p.7|) we reproduce the one-loop 
result which is a result of the perturbative unitarity of the chiral expansion. 
As a result when the sum rule for this combination is evaluated with the phys- 
ical absorptive parts, we get an answer that is substantially different from its 
one-loop expression. From this we conclude that the two-loop corrections to 
this combination must account for this discrepancy, although we cannot con- 
clude which one of the elements in the combination receives the correction. 
An analogous exercise may be performed for the 7r 7r° combinations and in 
particular for 30a 2 — b where a 2 = (a° + 2al)/3 and b = + 26q)/3. A 
variety of combinations that arise from totally symmetric amplitudes that 



have similar properties has been examined recently [27]. 

A final exercise we perform is to compute the value for 2a{j — 5<2q from the 
Roy equation fit and then to insert the value of U we find into the one-loop 
formula for the same combination: 

2<z°-5a2 = — „ n (h + 5/8)1- 

Once again we observe that the agreement, while being fair is not exact 



reflecting the somewhat large uncertainties in our determination referred to 
earlier as well as due to the 0(p 6 ) corrections to its one-loop formula. 

4 Numerical Results and Discussion 

The numerical solutions of the Roy equations are obtained by using a param- 
eterization similar to the one described in great detail in [18|. In this study 



we have employed the tttt scattering from the CERN-Munich experiment and 



documented by Ochs |28[ in the region 19 < s < 60, and the high precision 



K14 experimental data p{| for the phase shift difference 5° — 8{. We have de- 



vised some checks on our computation in the following manner: we have first 
of all determined the best fit to the parameters of the parameterization re- 
ferred to above by requiring it to simultaneously yield the best fit to the data 
as well as satisfy the Roy equations in the domain of their validity The same 
parameterization is considered to be valid in the domain 60 < s < 110 as 
well. Then we have employed the Ochs data in this region and have required 
a best fit to the data in 4 < s < 110 and that it satisfies the Roy equation in 
4 < s < 60. In practice the data above 60 do not influence the parameters of 
the fit significantly: this in turn renders our numerical results rather stable 
since many of the quantities we compute are dominated by the low-energy 
behavior. The details of our work will be documented elsewhere [[30 . 



Our solutions require as an input parameter only a® and numerically 
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search for those solutions that minimize the discrepancy with respect to the 
data. In Fig. 1 we present our phase shift fits for 5®— S\ obtained for a[j = 0.19, 
0.20, 0.21 and also for the central value a{] = 0.26. The final result of the 
procedure above is an explicit representation for the lowest partial waves 
as functions of the energy in terms of a few parameters that optimally fit 
the data up to A = 110 and verifying the Roy equations in their domain of 
validity. 

The numerical values of the chq, Jq, I = 0, 2 and j3\ may now be computed 
from the explicit phase shift representation provided by the Roy equation fit, 
and l\ and I2 are extracted from eq. (|2.20 ) and ( |2.21| ) (the effects arising 



from the neglect above the cutoff A are disregarded; an unrealistic absolute 
saturation of the integrands above A can lower li by 0.13 and increase I2 by 
the same amount), and U from the appropriate comparison of eq. fl2.16|) and 



eq.(|2~T9|). 



The Roy equation fits are also used to compute the combinations of S- and 



P-wave scattering lengths and effective ranges of eq.(3 1 2)-(3 : 5) and the high 
energy tail and higher wave contributions computed from the driving terms as 
expressed in them. For the right hand sides of the sum rules for the D- wave 
scattering lengths, the driving terms for the Roy equations for the D- waves is 
not available in the literature. The S- and P- wave contributions are explicitly 
accounted for up to the cut-off and their high energy tail is disregarded 
in our tabulation of the results (an absolute (unrealistic) saturation of the 
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integrands above A = 110 yields an error of +0.7 x 10~ 4 on a® and an error 
of +0.07 x 10~ 4 and —0.22 x 10~ 4 on af). The only higher wave contribution 
arises from the /2(1270) |3T| and we have estimated its contribution from 



two inequivalent sets of sum rules for these available in the literature [£6|, ^7 
which give almost identical results of 0.54 x 10~ 4 for a° and 0.38 x 10~ 4 for 
a\ and have been included in the final results. 

We tabulate our results in Table 1 and 2. In Table 1 for a given input of 
do we report the results of our fit for eft, &q, 6q, a\,b\, where the last four 
are obtained from the sum rules of the previous section and the values of li, 
I2 and U computed from the dispersive relations. In Table 2 (a), 2 (b) and 2 
(c) we give values for the 9 combinations of threshold parameters of interest 
computed from the sum rules and the Roy equation representation and also 
their one-loop values obtained by inserting the values of l\ and I2 from the 
dispersive analysis. 

An inspection of this table reveals that the values of the combinations for 
18a\ — (2<2q — 5<2q) in both columns agree better than the values in the two 
columns for 18b\. Such an agreement is also seen to be better for 02 than it 
is for bo. We conclude therefore that the determination of l\ and I2 from the 
dispersive framework is better consistent with the one-loop expressions for 
the scattering lengths than it is for the effective ranges. 

As a final check we have compared our results obtained from the Roy 
equation representation with those obtained from a simple analytic param- 
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eterization of the type proposed by Schenk |32|] for the lowest wave phase 
shifts: 

tantfOO = P (s) Li + [bi - - °° + (a^) 3 ]( g - 4)/4l J = 0, 2 

[ (s/ - 4J/4 J sj - s 

tan^C) = PW^M + ft! ~ ( ! 4 )/4 ]( g ~ 4 )/ 4 >^- 

This parameterization employs the postulated normal threshold behavior 
and the properties of threshold expansion of the real part of the partial waves 
and elastic unitarity and incorporates features of the tttt interaction such as 
the position of the p, s p = 30, and that the I — S- wave passes through 
7r/2, So = 36 (and an unphysical S2 = —16). The remainder of the required 
inputs are obtained from the first column in each of the Tables 1, 2 and 3 in 
addition to the values of a(j and that correspond to these and evaluated 
the «o,7o,/ = 0,2 and f3\ with A w 50 [K-K threshold) and solved for l\ 
and For instance for the entries of Table 2 (b) we find l\ = —1.72 and 
I2 = 5.2. Thus this simple parameterization confirms our numerical findings 
to within a surprising 5% accuracy. However when this parameterization is 
inserted back into the sum rules for the threshold parameters, the agreement 
is good only to about 15%. 

We now remark on contrasting our determinations of l\ and I2 with those 
mentioned in the introduction. Our numbers for l\ and 1% are comfortably 
accommodated in the range that was first obtained from D- wave scattering 
lengths in Ref. and more or less accommodated in the range obtained 
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from the K14 decays. Our value for l 2 is significantly lower; the somewhat 
larger values of the I = D- wave scattering length used in Ref. than 
the numbers we find appear to be the cause. In the present work we have 
presented a clear correlation between the input value of a{j and the D- wave 
scattering lengths which points towards a 10% smaller value for than the 
central value of 17 x 10 -4 employed there. Furthermore, we are working in 
a specific dispersive framework where the one-loop expressions for the effec- 
tive subtraction constants are altogether likely to suffer some higher order 
corrections, which we have not attempted to analyze in this work. Similarly 
the other determinations cited are also susceptible to such corrections and 
it has not been possible to explain quantitatively what the origins of the 
discrepancy are. 

5 Conclusions 

The coupling constants l\ and l 2 of the one loop chiral expansion has been 
accurately determined at 0(p 4 ) precision from a Roy equation analysis of 
the existing Tin scattering data with e (0.19,0.21) predicted by standard 
chiral perturbation theory. A suitable dispersive framework is used to effect a 
comparison between the one loop chiral representation and the Roy equation 
phase shift representation of the tttt amplitudes to obtain 

J x = -1.70 ± 0.15 and l 2 « 5.0. 
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The result is consistent with the bounds obtained on the combination l\ + 
2l 2 in Ref. Certain ambiguities in determinations of these at 0(p 4 ) 

from (combinations of) threshold parameters is discussed. The numerical 
consistency of one-loop results for those involving certain scattering lengths 
with the new values of l\ and I2 and their values computed from the Roy 
equation fits is superior to the consistency for those that involve effective 
ranges. 
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Note added 



After this work was completed we received an e-print [33] dealing with the 
subject of phase shift data, sum rules and chiral coupling constants at 0(p 6 ). 
The coupling constants in the first row of Table 4 therein corresponds to 
l\ = —1.36 and l 2 = 5.2. This may be compared with the coupling constants 
corresponding to our results for a® = 0.20. We find agreement to within a 
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few percent with the Roy equation solution while the agreement is very good 
for l 2 computed from Schenk's model. 



Appendix A 



Independent of the dynamics governing the interactions, it has been rigor- 
ously established that fixed-t dispersion relations with two subtractions may 
be written down for the amplitudes of definite isospin in terms of unknown 
t-dependent subtraction functions: 

T'(M,«) = j:C I J'(C I '(t) + (s-u)D I '(t)) + 



1 r°° dx 



7T J A 



-I 



II' 



U 



X* \ x — s 



x — u 



C 1 ^ A r {x,t), (A.l) 



where A\x, t) is the isospin / s-channel absorptive part, C st and C su are the 
crossing matrices: 



C. 



St 



/1/3 1 5/3 \ 

1/3 1/2 -5/6 
U/3 -1/2 1/6 ) 



/ 1/3 -1 5/3 \ 

-1/3 1/2 5/6 
V 1/3 1/2 1/6/ 



and I is the identity matrix. Bose symmetry implies: C 1 (t) = D°(t) = 
D 2 (i) = 0; the unknown t-dependent functions C / (t) and -D 7 (t) may be elim- 
inated, using crossing symmetry, in favor of the S-wave scattering lengths. 

The Roy equations are obtained upon projecting the resulting dispersion 
relation onto partial waves and inserting a partial wave expansion for the 
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absorptive part. They have been rigorously proved to be valid in the domain 
—4 < s < 60. These are a system of coupled integral equations for partial 
wave amplitudes of definite isospin / which are related through crossing sym- 
metry to the absorptive parts of all the partial waves. The Roy equations for 



the S- and P- waves are [16, 17, 18 



f° (s) = a° + (2al - 5a*)— - + E E / dxK 1 ^ (*, x)lmf} (x), 

e A ^ °° roc 

fl(s) = (2al-hal) S — + EE/ dxK['(8,x)1mf}{x), (A.2) 



2 oo 



e 4 " z~ roo 

f (s) = ag-(2a°-5ag)^ r +EE/ 4 dxK l J (s,x)lmff (x) 



s-4 

+ 

J'=OZ'=0' 

and for all the higher partial waves written as: 

#00 = EE/ Mj7'( S ,x)Im/^(x), />2, 

where Kj'/ (s, s') are the kernels of the integral equations and have been 
documented elsewhere [I7|. Upon cutting off the integral at a large scale A 
and absorbing the contribution of the high energy tail as well as that of all 
the higher waves over the entire energy range into the driving terms dj (s, A) 
we have: 
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f (s) = a° + (2a° - hal) S —^ + £ £ dxK™ (s, x)Tmff[x) 
+d? {s,A), 

fl(s) = (2a° -5al) S ^ + i:j: [\xK({(s,x)lmff 



X) 

I'=Ol'=0' " 

+d{{s,A), (A.3) 

_4 2 1 A 



/ 2 (s) = al-(2a° -5al) S —-+Y,T,f dxK^^Tmf}' 

Z4 /'=0/'=0' / ' 1 



X 



and for all the higher partial waves written as: 

fi(*) = EE/ A ^(^W^) + ft,A), ^>2. 

r=oz'=o 4 

The driving terms themselves for the two lowest partial waves, the I = 0, 2 
S- waves and the / = 1 P- wave are available in the literature when A = 110 
[15]. [We note that the driving terms for the I > 2 partial waves are not 
documented in the literature.] The primary aim of this work is to derive as 
completely as possible the information on the S- and P- waves and that of the 
numerical fit to the experimental data (described in a subsequent section) to 
provide a parametric representation for the physical S- and P- wave partial 
waves. 

Out of the absorptive parts of the physical S- and P- waves, one may 
construct a manifestly crossing symmetric amplitudes [O, 2j^: 
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32tt v ' ' ' 

J'=0 4 71 ^ 4 



<ix 



x — s 



■r r + 



X — t 



4 x(x — 4) 
C. 



x — « 



+3 



S(*-«) T 7/' , *(S-«)^,JI' , U(t-S 



-I 



x — s 



x — t 



x — u 



-c 



II' 



Im/o'(x) 
Imf('(x) 



(A.4) 



Our objective may be met by first rewriting this dispersion relation as: 



327T 



I'=0 



+(s(t - u)l ir + t{s - u)Ctf + u(t - s)C^)(3{ 



II' 



ill 1 \al' 



1 Ui a -r-r-rl 



7T J4 X J \ I - S 



x — u 



7r J4 x 2 (j; — 4) \ a; — s x — t 

where a^, Pi and 7q are defined in eq. (|2.18| ). We now invert the isospin 
amplitudes eq. fl2.2| ) to obtain A(s, t, it) that is constructed out of the S- and 
P- wave absorptive parts: 



A(s, t, u) 



1 



T°(s,t,u)-T 2 (s,t,u 



(A.5) 



and is given in eq. ( p.!7| ). 
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Table Captions 

Table 1. The computed values corresponding to the input a° of a^, b^, a\, 
b\ and the computed values of l\, I2 and U. 

Table 2. (a) Values of combinations of threshold parameters corresponding 
to the fit of the first line of Table 1 and their one loop values with the new 
li and l 2 ; (b) As in (a) for the second line; (c) As in (a) for the third line 

Figure Caption 

Fig. 1. Results of our fit to the combination 5$ — 5{ (in degrees) as a function 
of energy (MeV): full line a(] = 0.26 dash-dotted line: a° = 0.21, dashed line: 
a[] = 0.2 dash-double-dotted line: a° = 0.19. Also shown is the Rosselet 
data. 
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0° 

Oq 


a 2 Q 


6° 


6 2 


a\ 


b\ 


h 


h 


k 


1 


0.19 


-0.040 


0.238 


-0.074 


0.035 


0.006 


-1.80 


4.98 


0.87 


2 


0.20 


-0.037 


0.237 


-0.074 


0.035 


0.006 


-1.69 


4.97 


1.09 


3 


0.21 


-0.035 


0.238 


-0.075 


0.035 


0.006 


-1.58 


4.96 


1.30 



Table 1 





Roy equations 


One loop formula 


18a} - (2a° - 5ag) 


0.0459 


0.0472 


186} 


0.1079 


0.0788 


a% 


15 x 10" 4 


13 x 10" 4 


a 2 2 


0.6- 10~ 4 


0.4 • 10" 4 


36° - (2a° - 5ag) 


0.1326 


0.0924 


Qb 2 + (2a° - 5ag) 


0.1378 


0.1370 


a 2 


0.0005 


0.0005 


bo 


0.030 


0.026 


2a° - 5ag 


0.580 


0.561 



Table 2 (a) 
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Roy equations 


One loop formula 


18a} - (2a° - 5ag) 


0.0410 


0.0436 


186} 


0.1066 


0.0753 


a° 


15 x 10~ 4 


13 x 10~ 4 


a\ 


0.7- 10~ 4 


1.0 • 10" 4 


36° - (2a° - 5ag) 


0.1269 


0.0989 


66g + (2a° - 5ag) 


0.1423 


0.1400 


a 2 


0.0005 


0.0005 


6 


0.030 


0.027 


2a°, - 5ag 


0.585 


0.563 



Table 2 (b) 
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Roy equations 


One loop formula 


18a} - (2a° - 5ag) 


0.0360 


0.0424 


186} 


0.1053 


0.0740 


a° 


15 x 10~ 4 


13 x 10~ 4 


a\ 


0.7- 10~ 4 


1.0 • 10" 4 


36° - (2a° - 5ag) 


0.1200 


0.1008 


66g + (2a° - 5ag) 


0.1472 


0.1407 


a 2 


0.0005 


0.0005 


6 


0.030 


0.027 


2a°, - 5ag 


0.595 


0.567 



Table 2 (c) 
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